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The response of an interferometer changing its orientation with respect 
to a fixed reference frame is analyzed in terms of the beam-pattern factors 



' and the polarization-averaged antenna power pattern. Given the antenna's 

motion, the latter quantity describes the antenna's directionality as a function 
^ , of time. 

On \ An interesting case is represented by the class of motions where the de- 

' tector's plane is constrained to move on the surface of a cone of constant 

Tlj- ' aperture; at the same time, the two arms are rotating around a vertical axis. 

^ ■ This picture describes, in particular, the motion of LISA, a proposed space- 

, based laser interferometer, as well as of other planned missions. The overall 

\ sky's coverage, and that of the galactic plane in particular, is provided as a 

function of the cone's aperture. 

Similarly, one can consider the case of an earth-based interferometer. Us- 
ing the same formalism, one can derive a simple expression for the antenna 
pattern, averaged over the time of arrival of the signal, as a function of the 
position and orientation on the earth's surface. In particular, there turn out 
d \ to be two particular values for the terrestrial latitude and the inclination an- 

gle with respect to the local parallel which render the time-averaged antenna 
response perfectly isotropic. 

In the frequency domain, the general result is that the detector's motion 
introduces in the instrumental response to a long-duration continuous signal 
a few harmonics of the orbital frequency, whose magnitude depends on the 
source's position in the sky. In particular, we describe LISA's response to 
circularly polarized sinusoidal waves coming from a few known binary systems 
in our Galaxy. 
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I. INTRODUCTION 



Gravitational waves in the low-frequency regime (10~^ to 10~^ Hz) can only be observed 
from space, due to terrestrial disturbances. In space, the only technique currently avail- 
able, besides pulsar timing, is based on Doppler tracking of an interplanetary spacecraft 
However, this relatively inexpensive method has not provided enough sensitivity, thus 
far, for a detection. While better sensitivities may be expected in the near future, with 
advanced spacecraft such as CASSINI, much more ambitious projects for gravitational wave 
observatories in space have been proposed. Among these, the most promising detectors are 
based on space-born laser interferometry. In particular, LISA (Laser Interferometer Space 
Antenna) @J3, and OMEGA (Orbiting Medium Explorer for Gravitational Astrophysics) 0] 
consist of six drag-free, laser-bearing spacecraft, launched in orbit around the sun (LISA) or 
the earth (OMEGA). The six spacecraft would be placed, in pair, at the vertex of a triangle 
with 5 X 10^ km sides for LISA, and 5 times smaller for OMEGA. At each corner, the two 
spacecraft are phase locked through the exchange of a laser signal, replacing in this way the 
central mirror of an ordinary Michelson intereferometer. Each of the two probes sends a 
laser beam to a probe at each of the other two equilateral points, where the tracking signal is 
transponded back by phase-locked lasers, and the returning beams are eventually interfered. 

In order to keep the triangular constellation as stable as possible, elaborated orbits have 
been designed. In LISA case each spacecraft is orbiting a circle of radius 3 x 10^ km over a 
period of 1 yr. The interferometer plane, at an inclination of 60'' with respect to the ecliptic, 
is also rotating around the sun with the same periodicity. In OMEGA case the orbital plane 
is almost coincident with the ecliptic, and the interferometer is rotating around itself with 
a period of 53.21 days. 

The complicate motion is reflected in the time evolution of the interferometer's response 
to a source located in a fixed position in the sky. We will investigate the behavior of the 
antenna response in presence of a generic motion, and apply our results to the specific 
motions of interest. As a side-product of our analysis, we can also examine a terrestrial 
interferometer, where the motion is simply related to the earth's rotation around its axis, 
and study its antenna pattern as a function of the location and orientation on the earth's 
surface. 

In this section, we briefly recall the formalism describing the antenna response to a 
gravitational wave passing by, in the long- wavelength approximation P]. 
First, we introduce the wave symmetric trace-free (STF) tensor 



where the (complex) vector fh is defined in terms of the polarization vectors ex and eV as 



The tensor W represents the wave field as measured in the interferometer's proper rest 
frame. Then we define the STF detector tensor 
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where is the unit vector along the i-th arm. The interferometer response is the scalar 
obtained from the contraction of the wave tensor W with the detector tensor D 

R{t) = Wi,D'' = F+h+ + F^K. (4) 

The beam-pattern factors and Fx depend on the antenna's orientation with respect to 
the wave's propagation direction and polarization axes. 

We can choose the reference frame as in Fig. 1, with the x-axis of the {x,y,z) frame 
bisecting the interferometer's arms, so that the only non-null components of D in this 
reference frame are 



Di2 = D21 = sin(2n) , (5) 

where 2Q is the aperture angle. Therefore, in order to obtain R{t), we just need the com- 
ponent W12 in this particular frame. Fig. 1 also shows the Euler's angles 6, 0, and ip which 
transform from the interferometer's frame {x,y,z) to the wave reference frame {X,Y,Z). 
The latter is defined with the Z-axis opposite to the propagation direction, and the X and 
Y axes along ex and eV, respectively. 

It is easy to find, for the + polarization 

F+ = sin(211)[cos^cos(20) sin(2V^) 

+ ^ (1 + cos^ e) sin(20) cos(2V^)] . (6) 

The beam-pattern factor Fx is obtained from Eq. with the substitution ^ ip + tt/A, a 
well known polarization property of gravitational waves. However, when averages over the 
polarization angle ip are considered, we can assume, without loss of generality, = = h. 
The quantity of interest is thus the polarization-averaged antenna's power pattern 

which for the interferometer in Fig. 1 reads 

1 



P 



1 + cos^ ^ + 6 cos 6 — sin 9 cos(40) . (8) 



For future reference, note that this definition of the antenna pattern is not normalized to 
unity, the average of P{6,(f)) over the whole sky being 2/5. A plot of P{6,(f)) in polar 
coordinates is shown in Fig. 5a. 

Eq. (H) gives the instantaneous power pattern for a wave impinging from the direction 
{6, (p) in the interferometer's reference frame. If the detector is moving with respect to the 
source, then, apart from Doppler effects considerations, all we need to do is simply replace 6 
and (p with the appropriate functions 6{t) and (p{t). For example, if the antenna is rotating 
around its vertical axis with angular velocity u (Fig. 2), then we can obtain the antenna 
pattern at any time t from Eq. (§), with the substitution 

<p^<p-^o-u;t, (9) 
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where is some initial angle. As a matter of fact, this very simple case describes, with 
good approximation, the time evolution of OMEGA 0]. As mentioned in the Introduction, 
the idea of OMEGA is essentially similar to that of LISA, except that the six spacecraft are 
launched into a circular earth orbit, beyond the Moon orbit. The triangle has now 10^ km 
sides, and it is rotating around itself with a period of ~ 53 days. 

However, to mask the motion of the antenna with the apparent motion of the source is 
not always convenient, especially when dealing with a large number of sources, or when the 
motion is very complicate. In this paper, we shall introduce a more useful representation, 
where every quantity is referred to a fixed reference frame, so that the source's polar coor- 
dinates 6 and (f) remain constant, and the antenna response depends on time through the 
actual motion of the interferometer. 



II. ANTENNA PATTERN FOR A GENERIC MOTION 

We now introduce an arbitrary reference frame {x',y',z'), with the only requirement to 
be stationary with respect to fixed starsQ. To be more explicit, when dealing with a space- 
born interferometer, we can adopt an Ecliptic coordinate system. In the last section we will 
also consider a terrestrial interferometer, which is most easily described in an Equatorial 
frame. 

The full description of the antenna response requires six Euler's angles, defined as in 
Fig. 3. The orthogonal transformation from the wave's frame {X, Y, Z) to the fixed frame 
(x', y', z') is given by the orthogonal matrix 

(cos^cos?/^ — cos sin sin — (cos^sin-?/; + cos ^ sin cos t/^) 
sin (f) cos ijj + cos 9 cos sin — sin sin -0 + cos 9 cos cos ijj 
sin 9 sin -0 sin 9 cos ip 

The matrix B which transform from {x,y,z) to {x',y',z') is analogous to the matrix A, 
with the Euler's angles 9, 0, -0 replaced by the corresponding ones (, rj, ^. Thus, the complete 
transformation from the wave's frame to the detector's one is given by ■ A. Actually, as 
we may expect from Fig. 3, the angles and 1] appear in our results only in the combination 
S = (j) — rj. 

After a rather lengthy calculation, one ends up with the following expressions for the 
beam-pattern factors and 

F+ = sm{2Q) [Acos(20 cos(2V^) + B cos{2^) sin(2V^) + 

+ C sin(20 cos(2^/') + D sin(2^) sin(2^)] , (11a) 

Fx = sin(2fi) [5cos(20 cos(2^) - Acos{2^) sin(2V^) + 

+ D sin(20 cos(2V^) - C sm{2^) sin(2^)] . (lib) 

The coefficients A,B,C, and D in Eqs. (|TT|) depend only on the angles (, 9, and 6. They are 
explicitly given in Appendix A. 



sin 9 sin \ 
— sin 6* cos . (10) 
cos^^ / 



^Here and in the following, when we say 'stationary' (or 'fixed') we mean stationary over the 
characteristic time scale of the detector's motion. 
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We recall that the power pattern P is obtained squaring, and averaging over the polar- 
ization angle ip, the interferometer response. From Eqs. @, (0), and ( ITll) we eventually 
obtain 



P =J2{i^n + l^n cos(40] cos{n6) + a„sin(40 sm{n6)} , (12) 

n=0 

where now the coefficients Xn, fJ^n, and cr„ depend only on 9 and (. These coefficients, 
given in Appendix B, are quite complicate trigonometric polynomials of their arguments. 
Nonetheless, Eq. (|T2|) turns out to be very useful in practice. In fact, for the planned 
detectors considered in the present paper, the angle ( is constant, and thus the only possible 
time evolution is related to the sinusoidal functions of ^ and 6 which appear explicitly in Eq. 
(|l^). Before analyzing in more detail the proposed space-born interferometers, we consider 
a trivial application of Eq. (P^ ) . 



A simple example 



As a ffist test of Eq. (jT^) we can consider, as we did at the end of Sec. I, an interferometer 
which is rotating around its z-axis (see Fig. 2), so that its trivial motion is described by 



C = ^ = 0, 



^ = ^0 + UJt 



(13) 



From Eq. (0), and Eqs. (B.1)-(B.15) of Appendix B, we find that the only non-null coeffi- 
cients in Eq. ( |T^ ) are 



" 8 

/i4 = 0"4 



1 + cos^ 9 + 6 cos 9 
1. . „ 



-- sm 



(14) 
(15) 



so that 



P = - 



1 + cos^ 9 + 6 cos^ 9 - sin^ 9 cos [4(0 - - ^t)] 



(16) 



As pointed out in Sec. I, this result can be obtained much more easily directly from Eq. (^, 
with the substitution (^. Eq. (|l^) gives, with good approximation, the antenna pattern of 
OMEGA, with ~ 1.4 x 10"^ sec~\ 



III. SPACE INTERFEROMETERS: CONICAL MOTION 

In the previous section, we have considered the antenna power pattern associated to an 
unspecified motion of the detector. We will now focus our attention to the case of a space 
interferometer, which presents, independently of the particular project under investigation, 
some very general and interesting properties. 

Inserting the space interferometer in its orbit and keeping the interferometer configu- 
ration stable over the mission lifetime - at least two orbital periods - is a very demanding 
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navigation task, due principally to the perturbations from the earth and the other planets. 
For instance, one of the consequences of the complicate orbit is the fact that we cannot 
maintain equal distances between the probes. In a recent paper [Q, we modeled the noise 
that is introduced into the differenced data because of the unequal arms, and showed that 
the final accuracy of the interferometer is not compromised. Another example of the prob- 
lems we may face in a space-born interferometer is that, due to the earth disturbances, high 
Doppler rates would result. Hellings et al. |^ described a laser transmitter and receiver hard- 
ware system that provides the required readout accuracy and implements a self-correction 
procedure for the on-board frequency standard used for laser phase measurement. 

We will now consider the interferometer's orbit, and discuss its implication on the antenna 
response to a wave coming from a given direction in the sky. In particular, as we mentioned 
in the Introduction, LISA will orbit the sun at the earth's distance, as far behind the 
earth as possible. The plane containing the six probes, during its orbit, will remain always 
tangent to the surface of a cone of 60° aperture, and the detector itself will rotate in this 
plane with same periodicity - one year - but opposite direction. Fig. 4, reproduced from 0, 
shows LISA configuration. 

In this section, we will consider a LISA-like motion, characterized by a generic cone 
aperture. In other words, the motion of each of the three interferometers is assumed to be 
described by0 

C = const , (17a) 
r] = rio + ujt 5 = (f) - T]o - cut) , (17b) 

^ = ^o-ujt. (17c) 

Note that rj and ^ are counter-rotating. Note also that the three interferometers formed 
by the triangular configuration have initial values which differ from each other by 120°, 
whereas rjo is the same for all of them. For the sake of conciseness, since we are concerned 
here with the antenna power pattern, we will consider only one of the interferometers, leaving 
the possibility of exploiting the polarization sensitivity to future works. 

The proposed LISA orbit has ( = 60°, a critical value for the stability of the triangular 
configuration. Since the behavior of the antenna, in terms of sky's coverage, directionality, 
etc., is very sensitive to the inclination, we will keep C as a free parameter throughout this 
paper, and refer to Eqs. (|T7|) as describing a 'conical' motion. For a given (, the coefficients 
A„, fin, and cr^ are now functions of 9. They are explicitly given in Appendix C for the LISA 
case. 

In order to determine the sky's coverage during the detector's lifetime, we need to consider 
the time average of the antenna pattern P over one orbital period, which gives 

{P)t = Ao + i (/i4 + ^4) cos[4(0 -r]o- ^o)] (18) 



^we could eliminate one of the two initial conditions or Vo by simply rescaling the time, taking 
for instance the origin of time at the passage through the line of nodes (.^o = 0) or through the 
vernal equinox (770 = 0). 
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Fig. 5 shows a plot of {P)t in polar coordinates for various values of the cone's aperture C- 
Note that C = corresponds to an interferometer fixed in space; since we are taking ^ and 
7] counter-rotating, the interferometer does not rotate at all. 

We know from Eq. (j^) that an interferometer fixed in space can never detect waves im- 
pinging from four specific null-directions, given by 9 = 7r/2, and (p = kir /2 {k = 0, . . . ,3). 
We can now analyze what happens to these null-directions in the generic conical case, fo- 
cussing our attention to the ecliptic plane 6 = 7r/2. It is easy to show, from the behavior of 
the functions Aq and |/i4 -(- cr4|/2 (see Fig. 6), that there is a tendency for the null-directions 
on the ecliptic plane to remain visible, although the magnitude of this effect is strongly 
affected by the value of (: for C = the null-directions are obviously completely preserved, 
while for ( = 7t/2 the 0-dependence is very poor. LISA is much closer to the latter case, 
and actually its {P)t is almost independent on 9 as well, as we shall see in a while. 

Note that, if ^ and t] were co-rotating, instead of counter-rotating as in LISA, then in 
the last term of Eq. ([T^) we would have to make the substitution 

^ (/i4 + cr4) ^ ^ (/i4 - C^4) • (19) 

which means that the 0-term would be even smaller, in absolute value, compared to Aq. 

One might infer from Fig. 5d that, if the inclination is close to ( = n/2, then the zeros of 
the antenna pattern can be found in the direction orthogonal to the ecliptic plane. However, 
this is not the case, since ^ = implies that P is still given by Eq. (H), with 9 and replaced, 
respectively, by ( and ^. Thus, even for ( = 7^/2 the time average is non-zero at the poles. 

Another important issue related to Eq. ([ISD is the time-averaged antenna directionality. 
As we can predict from Fig. 5, directionality is strongly dependent on the angle (. To make 
this statement more precise, let us consider the r.m.s. deviation from isotropy, defined as 



A = 



175 
iSvr 

47r 



{{P)T-2/hf dn 



1/2 

(20) 



The normalization factor in front of Eq. (^) is chosen in such a way that A is normalized 
to one for a fixed interferometer. For a conical motion, inserting Eq. (plSf ) in Eq. (^OD we 
easily obtain A = A(C), shown in Fig. 7. Explicitly, we find 

A(C) = — [19779 + 120 cos C - 118380 cos^ C 
288 '- 

+840 cos^ C + 180690 cos^ C + 840 cos^ C 

-3180 cos^ C + 120 cos^ C + 2115 cos^ C] • (21) 

Thus, the antenna pattern is distributed more and more isotropically as we increase C, 
from C = After we reach a minimum at C — 55°, the antenna's directionality starts 
increasing again. 

Therefore, we can conclude that, in the conical case, it is impossible to get a perfectly 
isotropic response, i.e. we never get A = (see also the above discussion about the 0- 
dependence). However, we can get very close to this ideal situation, if we choose C appro- 
priately. Remarkably, LISA's inclination is very close to the optimal value C, — 55°, and 
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gives A(60°) ~ 0.14. This means that LISA, during its hfetime, will cover the whole sky in 
an approximately uniform manner. Of course, in some circumstances, directionality needs 
to be preserved. For instance, one may want to disentangle the isotropic component of the 
stochastic background from the anisotropic contribution of the galactic binaries. Direction- 
ality can always be preserved integrating over a shorter period. Fig. 8 shows LISA's antenna 
pattern averaged over 3 and 6 months, respectively. 

Finally, we want to consider the interferometer's responsiveness to the galactic plane, 
given that most of the strongest sources will lie on this plane. In our coordinate system, the 
galactic plane is characterized by 

9((j)) = arctan(a cos (f> + j3 sin , (22) 

where a ~ 1.75 and /5 ~ 4 x 10~^. Assuming, for simplicity, that the sources are distributed 
isotropically on the plane, the event rate of disk's sources is proportional to the average area 
of the intersection of this plane with the antenna pattern, equal to 



P{e{(j)),(j))fdtd(f), (23) 







^ 2n T 

where d{(j)) is the function given in Eq. (p2D. This area depends on the details of the detector's 
motion, in our case on (, ^o, and rjo. Fig. 9 shows the quantity G as a function of (, since 
the dependence on and rjo can be neglected in a first approximation. For LISA, we find 
the small value G ~ 0.46, the exact value depending on the initial conditions rjo + ^o- We 
conclude that LISA is not particularly sensititive to the galactic disk, due to the fairly large 
inclination to the Ecliptic of both the detector's plane and the Galaxy. According to this 
crude analysis, we can expect that OMEGA, with a smaller, almost negligible, inclination, 
would increase its chances of observing a signal from the disk by roughly a factor two. 
Note, however, that previous calculations have shown that gravitational radiation from 
galactic binaries in the disk is comparable to that coming from the local region (r < 200 
pc). Since the latter is isotropically distributed, the time- varying signal from the disk 
contributes only a fraction of the galactic binaries stochastic background. Nonetheless, this 
small component could make the stochastic signal distinguishable from the detector's noise. 
In a future paper we will investigate the amplitude modulations introduced by the antenna 
motion in the confusion noise generated by different populations of galactic binary systems, 
and describe how to exploit this effect in order to detect the signal and to obtain information 
about the distribution of sources in the Galaxy. 



IV. EARTH-BASED INTERFEROMETERS 

As an additional application of Eq. (|r^), let us consider a terrestrial interferometer. In 
this case the most convenient choice for the 'fixed' reference frame is the Equatorial one, 
with the z' axis directed toward the North Pole, and the x' axis toward the Vernal Equinox. 
In this frame, the motion of the interferometer becomes similar to the conical case previously 
analyzed, except that now the detector can not rotate on itself, of course. In other words 
we have 
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C = 7r/2-£, (24a) 
r] = rio + ujt, (24b) 
^ = t, (24c) 

where i is the terrestrial latitude, u is the earth's angular velocity of rotation, and l G 
[~f ' + l] is the angle between the arms' bisector and the local parallel. 

For example, let us consider the average of the antenna pattern P over the time of arrival 
of the signal. In the terrestrial CclSG, clS opposed to the conical case considered in Sec. Ill, 
the interferometer cannot rotate around its vertical axis, and therefore averaging over time 
or gives the same result, namely 

{P)t = Xo{0, i) + /io(^, i) cos(46) . (25) 

For any specific value of 6, for example 6 = 102°, corresponding to the direction of the 
center of the Virgo cluster, Eq. (|25D gives the square of the r.m.s. power as a function of 
the antenna's position and orientation on the earth's surface, a quantity already numerically 



studied in |jTO| 



Instead of fixing 9, we could try to answer the question: is there any particular location 
and orientation for which the antenna pattern, averaged over one day, is isotropic? The 
answer turns out to be affirmative, and the isotropic antenna is characterized by 

iis = arcsin (^±^) - ±35°.26438972 , (26a) 
Lis = ^ arccos (^-^^ ~ ±25".38423976 . (26b) 

One can easily check, by inspection, that iis and Lis produce {P)t = 2/5 or, equivalently, 
A = 0. A detector located at latitude iis and oriented by Lis maximizes the event rate of an 
isotropic population of sources. 

In the event that the detector's position has already been chosen, one can still make 



use of Eq. ([25| ) in order to find the optimal orientation t^, which gives the least directional 
antenna pattern at that latitude. At each latitude i, we define as optimal that orientation 
i=K which minimize the quantity A. Fig. 10 shows l^ = L^,{i) and the corresponding minimum 
A* = A(£, t*). 

We stress that, as one may actually expect, A depends much more strongly on i than 
L, and in particular the antenna becomes rapidly anisotropic as we move away from iis, 
no matter how optimally we try to choose l. Moreover, we can foresee several terrestrial 
interferometers to be operative in the near future, so that the sky's coverage of a single 
antenna is not really an issue as critical as in the space-born case previously discussed. 



V. SINUOSOIDAL WAVES FROM BINARIES. FOURIER ANALYSIS. 

It is generally assumed that galactic and extragalactic binary systems are the most 
promising sources of gravitational waves for detectors based on laser interferometry. In fact, 
waves from a binary star, including the effect of eccentricity, orbital inclination, and also 
post-Newtonian corrections, have long been studied, and are today well understood. In 
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particular, the sensitivity of the planned space interferometers should allow the detection of 
waves from several known galactic binary stars. In the LISA and OMEGA frequency band, 
the strongest among these sources are presumably the Interacting White Dwarfs Binaries 
(IWDB) Table I contains the available data for five IWDB, including the amplitude 

and frequency of the expected gravitational waves. We have applied our results to these 
objects, and describe the LISA's response to the waves originating from them. Fig. 11 
shows the beam-pattern factors and for the five IWDB in Table I, as seen from 
LISA over one year. Fig. 11 also shows the analogous quantities for a sinusoidal signal, of 
unspecified amplitude and frequency, coming from the galactic centre. 

We will now consider the effect of the motion in the frequency domain, for both the 
space-born and the terrestrial cases. We define the Fourier series as usual 

Pit) = «o + E cos 1^^^ + bk sin (^^^ I , (27) 

where T = 2L is the orbital period, and 

1 

ao = — jp{t)dt, (28) 

-L 

-L 



ak 



^ j P{t)coJ^\dt, (29) 



-L 

-L 



bk = 2 I Pit)^^^y-^jdt. (30) 

For the conical case analyzed in Sec. Ill, using eqs.(Bl)-(B14) of Appendix B, one finds that 
the only non-zero Fourier coefficients are given hj [Sq = (p — rjQ, k = 1, . . . ,8): 

ao = \o + ^(/^4 + cr4) cos[4((5o - ^o)] , (31) 

4 

Ofc = /ioCOs(4^o)'5fc4 + ^[XnCOs{nSo)Snk 



n=l 



+ ^(/i„ + cr„) cos{n6o - 4^o)^fc+n,4 



+ ^(/in - COs(n5o + 4^o)'5fc-n,4} , (32) 



bk = /iosin(4^o)44 + E{A„ sin(ra5o)5n^ 

n=l 

-^(/"n + (Tn) Sm{n6o - 4^o)4+n,4 



+^(yUn - o-„) sin(ra5o + 4^o)4-n,4} • (33) 



The analogous calculation for the terrestrial case, described by Eqs. (|2J), gives the 
following nine coefficients (A; = 1, . . . , 4) 
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oo = Ao + /io cos(4i) , (34) 
Ofc = [Afc + fik cos(46)] cos(Mo) + (Tk sin(4i) sm{kSo) , (35) 
bk = [Afc + /ifc cos(4i)] sin(A;5o) — o'k sin(4i) cos(fc5o) • (36) 

In general, given the initial detector's position, these Fourier coefficients depend on the 
source's coordinates 9 - through Xk,fik,o'k - and (p - through 6o. If the source location 
is known, then one can look for these spectral lines as a convincing signature about the 
gravitational origin of the signal. When the source's coordinates are unknown, however, one 



has to deal with the complication arising from the Doppler effect |]TT|. The motion of the 
detector, besides the amplitude modulation described in the present work, also introduces a 
location-dependent phase modulation, in the form of a Doppler broadening of the sinusoidal 
signal. In the case of LISA, the magnitude of this effect, over a period T = 1 yr, is 

^^2xl0-^ (37) 

so that, in the spectral region below 10~^ Hz, we do not have any hope of finding the 
aforementioned lines, separated from each other by only 1/T ~ 3 x 10~^ Hz. For a terrestrial 
interferometer, the situation is analogous, only complicated by the simultaneous effects of 
the diurnal and annual motion of the earth, and also by the earth-moon interaction. 

In conclusion, for long enough observations, we need special techniques to compensate 
for the frequency spread over several frequency-resolution bins, and eventually to recover the 
amplitude modulation described in this paper. Several different strategies can be adopted 
to overcome this problem, although none of them is completely satisfactory, due to the large 
amount of computation involved. See |1T[ for more details. In any case, the amplitude 



modulation can be exploited for an independent measurement of the location of the source, 
and, in addition, to obtain the polarization of the wave. 



VI. CONCLUSIONS 

In the present work, we have considered a gravitational wave interferometer, in motion 
with respect to fixed stars, and studied the resulting amplitude modulation of a long-duration 
continuous signal. The general results are presented in Sec. II, where the istantaneous beam- 
pattern factors - Eqs. (|TT|) - and the polarization-averaged antenna power pattern - Eq. (|T2|) 
- are given as functions of time, for a generic motion. 

Next, two particular cases have been analyzed: 1) the probable orbit of a space-born 
interferometer, with particular emphasis on LISA, and 2) the motion of a ground-based 
interferometer. 

For what concerns LISA, we have shown that its peculiar motion makes the time-averaged 
antenna pattern practically isotropic, thus providing an uniform coverage of the whole sky 
over the period of one year. For shorter integration periods directionality is mostly preserved, 
and can be exploited where necessary, for example in the search of a galactic binaries back- 
ground. However, when we focused on the galactic disk, we found that the average antenna 
response is far from optimal, due to the relative orientation of the Ecliptic and the galactic 
plane itself. We stress that these results are not conclusive, since we have neglected the 
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anisotropy in the distribution of the sources with respect to the sun, due to the fact that we 
are located near the edge of the disk. In this respect, additional work is needed. 

In the terrestrial case, thanks to the probable redundancy of future gravitational wave 
observatories, the discussion about a single antenna's sky's coverage is not so critical. How- 
ever, we found that there are particular positions on the earth's surface, given in Eqs. (p6[), 
which render the time- averaged antenna response perfectly isotropic. For what concerns the 
galactic plane, since the latter makes with the Equatorial plane approximately the same 
angle it makes with the Ecliptic (~ 60°), the result is analogous to the conical case, with 
the angle ( interpreted as 90 — £ in Fig. 9. In other words, the response to the galactic plane 
increases as we move the interferometer from the equator toward the poles, with a minor 
role played by the orientation angle l. 
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APPENDIX A: 

In this Appendix, we give the coefficients A,B,C, and D which enter in the expressions 
of F+ and Fx, Eqs. (|TT|). They, in turn, can be expressed in terms of intermediate quantities 
ai,a2,/3i, and j32 as follows 



A = ai/3i + 02/32 
B = ail32 - a2l3i 

D = aia2 + (3il32 



(Al) 
(A2) 

(A3) 
(A4) 



where 



ai = cos 6 

a2 = cos 9 sin S 

Pi = cos ( sin 6 

P2 = cos C cos 9 cos 5 + sin C sin 9 



APPENDIX B: 



In this Appendix, we give the coefficients Xn, fJ-n, and cr^ (n = 0, . . . , 4), defined in Eq. 
(|1^), as functions of the angles 9 and (. 
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Ao = ^ (35 + 35 cos^ e cos^ C - 30 cos^ ^ + 3 cos^ C + 108 cos^ C cos^ 9 

-30 cos^ C cos^ ^ - 30 cos^ C cos^ ^ - 30 cos^ C + 3 cos^ O) (Bl) 

Ai = ^ sm(2C) sin(2^) (l5 - 3 cos^ 9-2, cos^ C + 7 cos^ C cos^ 9) (B2) 

A2 = — (7 + 7 cos^ 9 cos^ C - 8 cos^ 9 + cos^ C + 16 cos^ C cos^ ^ - 8 cos^ C cos^ 9 

-8cos\cos^e-8cos^C + cos^^) (B3) 

A3 = - sin^ C cos C sin^ 9 cos 9 (B4) 
8 

A4 = -^sin^Csin^^ (B5) 

/io = - ^ sin^ C (3 - 30 cos^ ^ + 35 cos^ 9) (B6) 



/xi = sin(2^) cos C sin^ C (3 - 7 cos^ 9) (B7) 

//2 = ^ sin^ C sin^ ^ (l + cos^ C) (l - 7 cos^ 9) (B8) 

/i3 = sin(2C) sin^ 9 cos ^ (3 + cos^ () (B9) 

sin^ ^ (1 + cos^ C + 6 cos^ C) (BIO) 
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ao = (Bll) 

(7i = sin(2^) sin^ C (3 - 7 cos^ 9) (B12) 

(72 = ^ COS C sin^ C sin^ 9(l-l cos^ 9) (B13) 
8 ^ ' 

(T3 = - sin C cos 9 sin ^1 + 3 cos C) (B14) 
8 ^ ' 

(74 = -— sin^^cosC(l + cos2C) (B15) 



APPENDIX C: 

In this Appendix, we provide the coefficients A„, and (j„ (n = 0, . . . , 4) for LISA 

Ao = ^— f 443 - 37 cos^ ^ - 78 cos^ 9] (CI) 

" 1024 V ) ^ ' 

Ai = sin(2^) f57 - 5 cos^ 9) (C2) 
256 ^ ^ 

A2 = (9 - cos^ ^ - 8 cos^ 9) (C3) 

^3 = ^sin3^cos^ (C4) 
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A4 

fj,4 



9 . . 



1024 
9 

~1024 



sin^^ 



15 
256 



(3 + 35cos^^-30cos2 e) 
sm(2e) (3-7cos^ e) 



sin2^(l-7cos2 e) 



13^/3 . 3 



128 
41 



1024 



sin 9 cos 9 



sm^9 



(70 = 

(71 = sin(2^) (3 - 7 cos' 9) 

= ^^sin'e fl - 7 cos' 9) 



0-3 
(74 



64 

7v^ . ,3 



64 



sin 9 cos 6' 



128 



sin^ 9 



(C5) 

(C6) 

(C7) 
(C8) 

(C9) 
(CIO) 

(Cll) 

(C12) 

(C13) 

(C14) 
(C15) 
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FIGURES 



FIG. 1. The geometry of the interferometer. Note that the usual spherical polar coordinates of 
the source's position are (0,(/> — vr/2). 

FIG. 2. OMEGA: a space-born detector rotating around its vertical axis. Although the actual 
omega's motion is better described with the formalism introduced in Sec. Ill, this simple example 
can provide a good approximation to it. 

FIG. 3. The relation between the wave's {X,Y,Z), the detector's {x,y,z), and the fixed 
{x',y',z') reference frames. 

FIG. 4. LISA: a space-born interferometer in orbit around the sun. Reproduced from with 
permission. 

FIG. 5. The antenna pattern, averaged over time and polarization, for various values of the 
inclination angle C- a) C = 0. b) C = 30°, c) ^ = 60°, d) = 90°. LISA corresponds to case c). 

FIG. 6. The (^-dependence on the Ecliptic plane, as a function of the angle The bigger is 
the difference between Aq and |^4 + (T4I/2, the smaller is the dependence on the direction (p in the 
averaged antenna pattern. 

FIG. 7. The r.m.s. deviation from isotropy of the time-averaged antenna pattern, as a function 
of C. The minimum of A corresponds to the maximum attainable isotropic response. This minimum 
occurs at C =^ 55°, for which A ~ 0.08. For LISA, A(60°) ~ 0.14. 

FIG. 8. EISA's antenna pattern averaged over (a) 3 months and (b) 6 months. The initial 
conditions are Co = % = 0. 

FIG. 9. The response of a space interferometer to the galactic disk, as a function of the angle C- 
For simplicity, a small dependence on the detector's initial orientation, more precisely on rjo + ^Oj 
has been neglected. 

FIG. 10. The minimum r.m.s. deviation from isotropy A^, attainable at any given latitude. The 
corresponding optimal orientation is also shown. 

FIG. 11. The LISA normalized responses to the 5 IWDB systems shown in Table I, and to a 
hypothetical source in the galactic center, over one year of observations. The two curves give the 
beam-pattern factors -F+ (solid line) and Fx (dashed line). We have assumed the initial conditions 
Co = ^0 = 0, and a polarization angle ip = 2. Horizontal axis is time in years. 
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TABLES 



TABLE I. Data for 5 known IWDB. The first column is the name, the second and third column 
are, respectively, the angles 9 and (j) in the Ecliptic coordinate system. The last two columns gives, 
respectively, the predicted amplitude and frequency of the gravitational waves. 



Name 


e 

(degrees) 


(degrees) 


GW Amplitude 
(10-22) 


GW Frequency 
(10-3 Hz) 


AM CVn 


52.56 


260.38 


5.27 


1.94 


CR Boo 


72.10 


292.27 


2.82 


1.34 


V803 Cen 


120.31 


306.17 


0.89 


1.24 


CP Eri 


120.83 


151.77 


4.02 


1.16 


GP Com 


67.00 


277.73 


1.77 


0.72 
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